Nucleosynthesis: what direct reactions can do for it? 
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The reactions of relevance for stellar evolution are difficult to measure directly in the laboratory 
at the small astrophysical energies. In recent years indirect reaction methods have been developed 
and applied to extract low-energy astrophysical S-factors. These methods require a combination of 
new experimental techniques and theoretical efforts, which are the subject of this review. 



I. INTRODUCTION 

Fusion - Fusion cross sections can be calculated from 
the equation pQ 

a F (E) = 7rA 2 ^(2^+l)P,(S), (1) 

i 



where E is the center of mass energy, A = \Jh 2 /2mE is 
the reduced wavelength and i — 0, 1, 2, • • • . The cross 
section is proportional to 7rA 2 , the area of the quantum 
wave. Each part of the wave corresponds to different im- 
pact parameters having different probabilities for fusion. 
As the impact parameter increases, so does the angular 
momentum, hence the reason for the 2£+l term. Pe{E) 
is the probability that fusion occurs at a given impact pa- 
rameter, or angular momentum. Sometimes, for a better 
visualization, or for extrapolation to low energies, one 
uses the concept of astrophysical S-factor, redefining the 
cross section as 

a F (E) = ^S(E)exp[-27rr,(E)}, (2) 

where rj(E) — Z\Zie 2 jhv, with v being the relative ve- 
locity. The exponential function is an approximation to 
Po(E) for a square- well nuclear potential plus Coulomb 
potential, whereas the factor 1/E is proportional to the 
area appearing in Eq. [T] 

Many reaction channels - Eq. [I] does not work in 
most situations. Only by including coupling to other 
channels, the fusion cross sections can be reproduced. In 
coupled channels schemes one expands the total wave- 
function for the system as 

^ = ^2a l (a)cj) i {a,q k ), (3) 

i 

where <fi form the channel basis, a is a dynamical vari- 
able (e.g., the distance between the nuclei), and q k are 
intrinsic coordinates. Inserting this expansion in the 
Schrodinger equation yields a set of CC equations in the 
form [2] 

|*=£a*<&|Er|fc>e«W* (4) 
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where U is whatever potential couples the channels k 
and j and E a — E&' — E^ is some sort of transition 
energy, or transition momentum. In the presence of con- 
tinuum states, continuum-continuum coupling (relevant 
for breakup channels) can be included by discretizing the 
continuum. This goes by the name of Continuum Dis- 
cretized Coupled- Channels (CDCC) calculations. There 
are several variations of CC equations, e.g., a set of dif- 
ferential equations for the wavef unctions, instead of us- 
ing basis amplitudes. Coupled channels calculations with 
a large number of channels in continuum couplings are 
somewhat challenging: the phases of matrix elements can 
add destructively or constructively, depending on the sys- 
tem and on the nuclear model. Such suppressions or en- 
hancements are difficult to interpret. 

Radiative capture - For reactions involving light nu- 
clei, only a few channels are of relevance. In this case, 
a real potential is enough for the treatment of fusion. 
For example, radiative capture cross sections of the type 
n + x — s- a + 7 and ttL (tt — E, (M) =electric (magnetic) 
L-pole) transitions can be calculated from [5] 

oizU = const - x Wcic \\GiriW l b j b )\ 2 , (5) 

where O^l is an EM operator, and (l c j c \\0„l\\ Ibjb) is a 
multipole matrix element involving bound (b) and contin- 
uum (c) wavefunctons. For electric multipole transitions 
(0^ L =r L Y LM ), 

POO 

{lcjc \\Oel\\ hjb) = const, x dr r L u b (r)u c (r), (6) 

Jo 

where Ui are radial wavefunctions. The total direct cap- 
ture cross section is obtained by adding all multipolarities 
and final spins of the bound state (E = E nx ), 

o- d -<-(E) = J2( SF ^» a L%W> ( ? ) 

L,J b 

where (SF)j b are spectroscopic factors. 

Asymptotic normalization coefficients - In a mi- 
croscopic approach, instead of single-particle wavefunc- 
tions one often makes use of overlap integrals, Ib{r), and 
a many-body wavefunction for the relative motion, u c (r). 
Both I b (r) and u c (r) might be very complicated to calcu- 
late, depending on how elaborated the microscopic model 
is. The variable r is the relative coordinate between the 
nucleon and the nucleus x, with all the intrinsic coor- 
dinates of the nucleons in x being integrated out. The 
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direct capture cross sections are obtained from the cal- 



culation of o^- j b oc 



\(h(r)\\r L Y L \\* c (r)) 



The imprints of many-body effects will eventually dis- 
appear at large distances between the nucleon and the nu- 
cleus. One thus expects that the overlap function asymp- 
totically matches (r — > oo), 



h{r) = Q X 



r 



(2nr) 



Ki b+ i/ 2 (Kr) 



,(8) 



where () are for protons and [] for neutrons. The binding 
energy of the n + x system is related to k by means of 
Eb = h 2 K 2 /2m n x, Wp^ q is the Whittaker function and 
is the modified Bessel function. In Eq. [HI Ci is the 
Asymptotic Normalization Coefficient (ANCjT^ 

In the calculation of cr^ °j above, one often meets the 
situation in which only the asymptotic part of h(r) and 
v I / c (r) contributes significantly to the integral over r. In 
these situations, u c (r) is also well described by a simple 
two-body scattering wave (e.g. Coulomb waves). There- 
fore the radial integration in crjj c j can be done accurately 
and the only remaining information from the many-body 
physics at short-distances is contained in the asymptotic 
normalization coefficient Cj, i.e. a^'j. oc Cf. We thus 
run into an effective theory for radiative capture cross sec- 
tions, in which the constants Ci carry all the information 
about the short-distance physics, where the many-body 
aspects are relevant [5j . It is worthwhile to mention 
that these arguments are reasonable for proton capture 
at very low energies, because of the Coulomb barrier. 

Resonating group method - One immediate goal 
can be achieved in the coming years by using the Res- 
onating Group Method (RGM) or the Generator Coordi- 
nate Method (GCM). These are a set of coupled integro- 
differential equations of the form [5] 



E 



J dV[<«(r,r')-M^(r,r')] to (r') = 0, (9) 



where H{N)(r,r') = (* A (a,r)| J ff(l)|* B (a',r')>. In 
these equations H is the Hamiltonian for the system of 
two nuclei (A and B) with the energy E, ^a,b is the 
wavefunction of nucleus A (and B), and g a (r) is a func- 
tion to be found by numerical solution of Eq. [9j which 
describes the relative motion of A and B in channel a. 
Full antisymmetrization between nucleons of A and B 
are implicit. Modern nuclear shell- model calculations are 
able to provide the wavefunctions ^ a b for light nuclei 

The astrophysical S-factor for the reaction 7 Be(p, 7) 8 B 
was calculated [7] and excellent agreement was found 
with the experimental data in both direct and indirect 
measurements [71 [8]. The low- and high-energy slopes 
of the S-factor obtained with a many-body microscopic 
calculation [7] is well described by the fit 



Sn(E) = (22.109 eV.b) 
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FIG. 1: World data on 7 Be(p,7) 8 B compared to theoretical 
calculations. 



where E is the relative energy (in MeV) of p+ 7 Be in their 
center-of-mass. This equation corresponds to a Pade ap- 
proximant of the S-factor. A subthreshold pole due to 
the binding energy of 8 B is responsible for the denomi- 
nator [9l[T0]. Figure[l]show the world data on 7 Be(p,7) 8 B 
compared to a few of the theoretical calculations. The 
recent compilation published in Ref. |11| recommends 
Sn = 20.8 ± 0.7 (expt) ± 1.4 (theor) eV b. 



II. DIRECT REACTIONS AND THE ROLE OF 
RADIOACTIVE BEAMS 

Transfer reactions - Transfer reactions A(a,b)B are 
effective when a momentum matching exists between the 
transferred particle and the internal particles in the nu- 
cleus. Thus, beam energies should be in the range of a 
few 10 MeV per nucleon. Low energy reactions of astro- 
physical interest can be extracted directly from breakup 
reactions A+a — > b+c+B by means of the Trojan Horse 
method (THM) [12 . If the Fermi momentum of the par- 
ticle x inside a = (b + x) compensates for the initial pro- 
jectile velocity v a , the low energy reaction A + x = B + c 
is induced at very low (even vanishing) relative energy 
between A and x. Basically, this technique extends the 
method of transfer reactions to continuum states. Very 
successful results using this technique have been reported 

Another transfer method, coined as ANC technique [H 
IT51 HE] relies on fact that the amplitude for the radiative 
capture cross section b + x — > a + 7 is given by 

M =< 4 a x (r bx )|0(r bx )|^ (+) (r bx ) >, 
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where 

i£ x =< <M&, £*, i-bx) > 

is the integration over the internal coordinates and 
^a:, of 6 and x, respectively. For low energies, the 
overlap integral I? is dominated by contributions from 
large r bx . Thus, what matters for the calculation of 
the matrix element M is the asymptotic value of 7™ ~ 
c bx W-r, a ,i/2(2K bx r bx )/r bx , where is the ANC and W 
is the Whittaker function. This coefficient is the prod- 
uct of the spectroscopic factor and a normalization con- 
stant which depends on the details of the wave function 
in the interior part of the potential. Thus, C£ x is the 
only unknown factor needed to calculate the direct cap- 
ture cross section. These normalization coefficients can 
be found from: 1) analysis of classical nuclear reactions 
such as elastic scattering [by extrapolation of the exper- 
imental scattering phase shifts to the bound state pole 
in the energy plane], or 2) peripheral transfer reactions 
whose amplitudes contain the same overlap function as 
the amplitude of the corresponding astrophysical radia- 
tive capture cross section. One of the many advantages 
of using transfer reaction techniques over direct measure- 
ments is to avoid the treatment of the electron screening 
problem |13[ I16j . 

Intermediate energy Coulomb excitation - The 
Coulomb excitation cross section is given by 

da.^f fdo\ 167r 2 Z|e 2 

dn ~ \dn) cl T? 

x E^tTF^I^m)! 2 , (ID 

7rXfi 

where B(irX,Ii — » If) is the reduced transition proba- 
bility of the projectile nucleus, irX = El, E2, Ml, ... is 
the multipolarity of the excitation, and /i = — A, — A + 
1 A. 

The relativistic corrections to the Rutherford formula 
for (dcr/dfi) , (relevant for collisions at 50 MeV/nucleon 
and above) has been investigated in Ref. [IT]. It was 
shown that the scattering angle increases by up to 6% 
when relativistic corrections are included in nuclear col- 
lisions at 100 MeV/nucleon. The effect on the elastic 
scattering cross section is even more drastic: up to 13% 
for center-of-mass scattering angles around 0-4 degrees. 

The orbital integrals S(ttX, /i) contain the information 
about relativistic corrections. Inclusion of absorption ef- 
fects in S(irX, h) due to the imaginary part of an opti- 
cal nucleus-nucleus potential where worked out in Ref. 
[18] . These orbital integrals depend on the Lorentz fac- 
tor 7 = (1 — ti 2 /c 2 ) -1 / 2 , with c being the speed of light, 
on the multipolarity 7rA/i, and on the adiabaticity param- 
eter £(&) = uifib/~fv < 1, where tofi — (Ef — E{) /fi is 
the excitation energy (in units of ft) and b is the impact 
parameter. 

Ref. [Tj5] has shown that at 10 MeV/nucleon the rel- 
ativistic corrections are important only at the level of 



1%. At 500 MeV/nucleon, the correct treatment of the 
recoil corrections is relevant on the level of 1%. Thus the 
non-relativistic treatment of Coulomb excitation [50] can 
be safely used for energies below about 10 MeV/nucleon 
and the relativistic treatment with a straight-line tra- 
jectory [21] is adequate above about 500 MeV/nucleon. 
However at energies around 50 to 100 MeV/nucleon, ac- 
celerator energies common to most radioactive beam fa- 
cilities, it is very important to use a correct treatment 
of recoil and relativistic effects, both kinematically and 
dynamically. At these energies, the corrections can add 
up to 50%. These effects were also shown in Ref. [55] for 
the case of excitation of giant resonances in collisions at 
intermediate energies. 

A reliable extraction of useful nuclear properties, like 
the electromagnetic response (B(E2)-values, 7-ray angu- 
lar distribution, etc.) from Coulomb excitation experi- 
ments at intermediate energies requires a proper treat- 
ment of special relativity [TU [23]. The dynamical rel- 
ativistic effects have often been neglected in the analy- 
sis of experiments elsewhere (see, e.g. [24]). The effect 
is highly non-linear, i.e. a 10% increase in the veloc- 
ity might lead to a 50% increase (or decrease) of certain 
physical observables. A general review of the importance 
of the relativistic dynamical effects in intermediate en- 
ergy collisions has been presented in Ref. [5J [55] . 

The Coulomb dissociation method - The Coulomb 
dissociation method is quite simple. The (differential, 
or angle integrated) Coulomb breakup cross section for 
a + A — > b + c + A follows from Eq. [Tl] It can be 
rewritten as 

da ^=N^(u J ;e-A).a;l a ^ b+c (u), (12) 

where u) is the energy transferred from the relative mo- 
tion to the breakup, and _y b+ c ( w ) is the photo nu- 
clear cross section for the multipolarity 7rA and photon 
energy u>. The function iV 7rA , sometimes called virtual 
photon numbers, depends on uj, the relative motion en- 
ergy, nuclear charges and radii, and the scattering angle 
Q = (0,cj)). N 1tX can be reliably calculated [26] for each 
multipolarity nX. Time reversal allows one to deduce 
the radiative capture cross section b + c — > a + 7 from 
CT 7+a -> 6+c( CJ )- This method was proposed in Ref. [57] 
and has been tested successfully in a number of reactions 
of interest for astrophysics. The most celebrated case is 
the reaction 7 Be(p, 7) 8 B [55], followed by numerous ex- 
periments in the last decade (see e.g. Ref. |29j). 

Eq. [12] is based on first-order perturbation theory. It 
also assumes that the nuclear contribution to the breakup 
is small, or that it can be separated under certain ex- 
perimental conditions. The contribution of the nuclear 
breakup has been examined by several authors (see, e.g. 
[50]V 8 B has a small proton separation energy (sa 140 
keV) . For such loosely-bound systems it had been shown 
that multiple-step, or higher-order effects, are impor- 
tant [31] . These effects occur by means of continuum- 
continuum transitions. Detailed studies of dynamic con- 
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tributions to the breakup were explored in refs. [351 [33J 
and in several other publications which followed. The role 
of higher multipolarities (e.g., E2 contributions [31 1341135] 
in the reaction 7 Be(p, 7) 8 B) and the coupling to high- 
lying states has also to be investigated carefully. It has 
also been shown that the influence of giant resonance 
states is small [35] , 

Knock-out reactions - The early interest in knock- 
out reactions came from studies of nuclear halo states, for 
which the narrow momentum distributions of the core 
fragments in a qualitative way revealed the large spa- 
tial extension of the halo wave function. It was shown 
[37j that the longitudinal component of the momentum 
(taken along the beam or z direction) gave the most 
accurate information on the intrinsic properties of the 
halo and that it was insensitive to details of the col- 
lision and the size of the target. In contrast to this, 
the transverse distributions of the core are significantly 
broadened by diffractive effects and by Coulomb scatter- 
ing. For experiments that observe the nucleon produced 
in elastic breakup, the transverse momentum is entirely 
dominated by diffractive effects, as illustrated 38 by the 
angular distribution of the neutrons from the reaction 
9 Be( 11 Be, 10 Be+n)X. In this case, the width of the trans- 
verse momentum distribution reflects essentially the size 
of the target [35] , 
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FIG. 2: Total knockout cross sections 
for removing the I — halo neutron of 
15 C, bound by 1.218 MeV, in the reaction 
9 Be( 15 C, 14 C 9S ). The solid curve is ob- 
tained with the use of free nucleon-nucleon 
cross sections. The dashed curve includes 
the geometrical effects of Pauli blocking. 
The dashed-dotted curve is the result us- 
ing the Brueckner theory, and the dotted 
curve is a phenomenological parametriza- 
tion. 

To test the influence of the medium effects in nucleon 
knockout reactions, we consider the removal of the I = 
halo neutron of 15 C, bound by 1.218 MeV. The reaction 
studied is 9 Be( 15 C, 14 C gs ). The total cross sections as a 
function of the bombarding energy are shown in figure [2] 
The solid curve is obtained with the use of free nucleon- 
nucleon cross sections. The dashed curve includes the 



geometrical effects of Pauli blocking. The dashed-dotted 
curve is the result using the Brueckner theory, and the 
dotted curve is the phenomenological parametrization of 
the free cross section. 
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FIG. 3: Longitudinal momentum distribu- 
tion for the residue in the 9 Be( 11 Be, 10 Be), 
reaction at 250 MeV/nucleon. The dashed 
curve is the cross section calculated using 
the NN cross section from the Brueckner 
theory and the solid curve is obtained the 
free cross section. 

In figure [3j we plot the longitudinal momentum 
distributions for the reaction 9 Be( 11 Be, 10 Be), at 250 
MeV/nucleon [40]. The dashed curve is the cross section 
calculated using the NN cross section from the Brueck- 
ner theory and the solid curve is obtained the free cross 
section. One sees that the momentum distributions are 
reduced by 10%, about the same as the total cross sec- 
tions, but the shape remains basically unaltered. If one 
rescales the dashed curve to match the solid one, the dif- 
ferences in the width are not visible [41;:. 

III. CONCLUSIONS 

There were many questions not addressed in this re- 
view, such as the role of central nucleus-nucleus colli- 
sions in determining phase transition, equation of state, 
and a quark-gluon plasma, all topics or relevance in as- 
trophysics. The review was more focused on the role 
of short-lived, exotic nuclei. There are many important 
scientific questions to be addressed both experimentally 
and theoretically in nuclear physics of exotic nuclei with 
relevance for astrophysics. These questions provide ex- 
treme challenges for experiments and theory. On the 
experimental side, producing the beams of radioactive 
nuclei needed to address the scientific questions has been 
an enormous challenge. Pioneering experiments have es- 
tablished the techniques and present-generation facilities 
have produced first exciting science results, but the field 
is still at the beginning of an era of discovery and ex- 
ploration that will be fully underway once the range of 
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next- generation facilities becomes operational. The the- 
oretical challenges relate to wide variations in nuclear 
composition and rearrangements of the bound and con- 
tinuum structure, sometimes involving near-degeneracy 
of the bound and continuum states. The extraction of re- 
liable information from experiments requires a solid un- 
derstanding of the reaction process, in addition to the 
structure of the nucleus. In astrophysics, new observa- 



tions, for example the expected onset of data on stellar 
abundances, will require rare-isotope science for their in- 
terpretation. 
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